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We study the nature of the inner Cauchy horizon of a Reissner-Nordstrom black hole in a quantum 
context by means of the horizon wave-function obtained from modelling the electrically charged 
source as a Gaussian wave-function. Our main finding is that there are significant ranges for the 
black hole mass (around the Planck scale) and specific charge for which the probability of realising 
the inner horizon is negligible. This result suggests that any semiclassical instability one expects 
near the inner horizon may not occur in quantum black holes. 
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I. INTRODUCTION 

Current literature contains many attempts at quantiz¬ 
ing black hole (BH) metrics, which focus on the purely 
gravitational degrees of freedom, and yield a description 
of the horizon unrelated to the matter state that sources 
the geometry [1]. However, this point of view might miss 
important features emerging from the highly non-linear 
nature of the gravitational interaction. One can think 
in analogy to the hydrogen atom, and note that its en¬ 
ergy levels, which have no classical counterpart, cannot 
be simply obtained by quantising the free electromag¬ 
netic field. A different perspective is taken in the Horizon 
Wave Function (HWF) formalism [2, 3], which is instead 
based on the quantum version of the Einstein equation 
relating the size of the gravitational radius (which can 
be a horizon) to the (quantum) state of matter. This 
formalism has been applied to several case studies [4, 5], 
yielding sensible results in agreement with (semi)classical 
expectations, and there is therefore hope that it will help 
our understanding of the quantum nature of BHs. 

In practical terms, the construction of the HWF starts 
from the spectral decomposition of the quantum mechan¬ 
ical state that represents a matter source localized in 
space. By expressing the energy in terms of the gravita¬ 
tional (Schwarzschild) radius, as it would be classically 
determined according to the Einstein equations, the spec¬ 
tral decomposition then directly yields the HWF. The 
normalised HWF supplies the probability for an observer 
to detect a gravitational radius of a certain areal radius, 
centred around the source in the quantum state that was 
used in the first place. The gravitational radius can then 
be interpreted as a horizon if the probability of finding 
the particle inside it, is reasonably high. In other words, 
the horizon size is necessarily “fuzzy” in this QM descrip¬ 
tion, just like the position of the particles that sources 
the geometry. 
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One would expect that quantum modifications of the 
sort mentioned above should be important only for small 
BHs, with masses close to the Planck scale, and that the 
quantum corrections to the classical solutions of large 
mass BHs should be negligible. A large BH can have 
arbitrary small curvature near its horizon, and we have 
numerous tests of gravity in such low density regimes. 
However, BHs are very peculiar objects, in that they trap 
any signal inside the horizon by definition, no matter how 
weak tidal forces are in its neighbourhood, and it may 
happen that they are better described as macroscopic 
quantum states (see for instance, the Bose-Einstein Con¬ 
densate model of Refs. [6]). 

Charged BHs were subject to many theoretical stud¬ 
ies in the past [7]. In this work, we shall extend the 
HWF formalism to the Reissner-Nordstrom (RN) geom¬ 
etry and investigate the probability for an electrically 
charged source represented by a Gaussian wave-packet 
to be a BH, and to have an actual inner horizon. In the 
semiclassical approach, the latter is a Cauchy horizon 
and is sometimes associated with an instability known 
as “mass inflation”: any small matter perturbation will 
blue-shift unboundedly just outside this horizon, and in¬ 
evitably produce a large deformation to the background 
geometry [8]. Although the existence of this effect is still 
debated (see, e.g. Refs. [9]), it is clear from the classi¬ 
cal causal structure of the RN geometry that, if there is 
matter falling through the outer horizon, it should ac¬ 
cumulate outside the inner Cauchy horizon, and even¬ 
tually lead to a large backreaction there. On the other 
hand, matter inside the inner horizon could escape the 
Cauchy horizon and produce a deformation as well. It 
is thus interesting to study under which conditions the 
inner horizon survives in the QM treatment. 

II. THE HWF FORMALISM 

The formalism introduced in Refs. [2, 3] is based on 
lifting the gravitational radius Ru of a QM system to 
the rank of a quantum operator. The coordinate r in a 
spherical metric is invariantly related to the geometrical 
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area At: r 2 of the sphere centred on the origin r = 0, 
therefore being a suitable candidate for an observable 
in the quantum theory. Moreover, the horizon radius 
r = i?H represents the location of trapping surfaces (sur¬ 
faces where the escape velocity equals the speed of light) 
and thus determines the causal structure of the space- 
time, which again is likely an observable property in the 
quantum theory. 

We remind the readers that in a neutral spherically 
symmetric system, 

f? H (r) = 2£ p ^M, (1) 

TO p 

where £ p is the Planck length and m v the Planck mass *, 
and 

M(r,t) = At: f p(f,t)r 2 dr (2) 

Jo 

is the Misner-Sharp mass. Here M = M(r,t) represents 
the total energy inside the sphere of area Anr 2 (thus, 
roughly speaking, including the negative gravitational 
energy) and is related to the energy density p of the 
source via the flat space volume. A specific value of r 
is a trapping surface if Rn(r) = r, while, if Rn(r) < r, 
the gravitational radius is still well-defined but does not 
correspond to any causal surface. 

Going back to the HWF formalism, let us start from 
QM states representing spherically symmetric objects, 
which are localized in space and at rest in the chosen ref¬ 
erence frame. Such particles are consequently described 
by wave-functions f/>s € L 2 (i? 3 ), which we assume can be 
decomposed into energy eigenstates, 

\^)=Y j C{ e )\^e) , ( 3 ) 

E 

where the sum represents the spectral decomposition in 
Hamiltonian eigenmodes, 

H I ) = E I 1>E) ■ (4) 

The actual Hamiltonian H needs not be specified yet 2 . 

The expression of the Schwarzschilcl radius in Eq. (1) 
can be inverted to obtain 

M = TOp It ’ ^ 

and remembering that M represents the total energy, it 
can be used to define the (unnormalized) “horizon wave- 
function” as 

V’h(-Rh) = C ( m p 4 ?h/ 2^ p ) • (6) 


1 We shall use units with c = kg = 1, and always display the 
Newton constant Cfv = t p /m p , so that h = t p m p . 

2 This is where, for instance, the self-gravity of the particle may 

enter. 


The normalisation of ipu is finally fixed by means of 
the scalar product 

/♦OO 

( V’H | </>H ) = 4 7T / ^(tf H )<M#H)i%di?H , (7) 
Jo 

again in agreement with the geometrical meaning of the 
variable as yielding the area 4 n R ^ of the correspond¬ 
ing sphere. Our interpretation of the normalised i/>h is 
then that it yields the probability for an observer to de¬ 
tect a horizon (necessarily “fuzzy”, like the position of 
the particle itself) of areal radius r = 7?h, associated 
with the particle in the quantum state ips- 

In more details, starting from the he wave-function tI>h 
associated with tl> g, we can now calculate the probability 
density for the particle to lie inside its own horizon of 
radius r = Ru- 

V < (r<R K ) = P s (r<R H )V ii (Rn), (8) 

where 

pRh 

P s (r < i? H ) = 4 tt / |V’s( r )| 2 f 2 dr (9) 

Jo 

is the probability that the particle is inside a sphere of 
radius r = i?n, and 

Ph(Ph) = 4 7t IV'h(Ph)| 2 (10) 

is the probability density that the sphere of radius r = 
i?H is a horizon. Finally, the probability that the particle 
is a black hole will be obtained by integrating (24) over 
all possible values of the horizon radius, namely 

/♦OO 

Pbh= V<{r <i? H )di? H ■ (11) 

Jo 

To anticipate the calculations in the next section, we 
note that electrically charged BHs have two gravitational 
radii, and we will therefore define two corresponding op¬ 
erators, namely R±. We shall then define a HWF for 
each one, and obtain the probabilities for both the inner 
and outer horizon to exist. 


III. ELECTRICALLY CHARGED SPHERICAL 
SOURCES 


We start from the RN metric, which can be written as 


ds 2 


with 


-/ d t 2 + r 1 dr 2 + r 2 (d6 2 + sin 2 6 dcj) 2 ) ,(12) 


2l p M + 
to p r r 2 


(13) 


where M and Q respectively represent the ADM mass 
and charge of the source. It is now convenient to intro¬ 
duce the specific charge 


\Q\m p 

£ P M 


( 14 ) 
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The case a = 0 reduces to the neutral Sclrwarzschild met¬ 
ric, and shall not be reconsidered here (see Refs. [2-5]). 
For 0 < a < 1, the above metric contains two horizons, 
namely 


R± = i p — ± 


=i P -( i±yr^ 
\ 



a 


(15) 


and represents a BH. The two horizons overlap for a = 1, 
the so-called extremal BH case, while for a > 1 no hori¬ 
zon exists and the central singularity is therefore accessi¬ 
ble to an outer observer. This is the prototype of a naked 
singularity, but we shall not consider this case here (the 
corresponding HWF is the topic of Ref. [10]). 

We shall now investigate the case 0 < a < 1, which 
classically possesses at least one horizon, from a QM per¬ 
spective. We first determine the HWFs and then calcu¬ 
late the probabilities for both the inner and outer horizon 
to exist. For this purpose, the classical relations (15) will 
be lifted to the rank of equations for the operators R± 
and M, which are chosen to act multiplicatively on the 
HWF (with the specific charge a viewed as a simple pa¬ 
rameter) . 


A. HWF for Gaussian source 

The source for the RN space-time is taken to be an 
electrically charged massive particle at rest in the ori¬ 
gin of the reference frame, represented by the spherically 
symmetric Gaussian wave-function 

r 2 

Mr) = p/2^/A ■ ( 16 ) 

The width of the Gaussian t is assumed to be the min¬ 
imum compatible with the Heisenberg uncertainty prin¬ 
ciple, that is 

e = x m ~ e p ^ , ( 17 ) 

TO 

where \ m is the Compton length of the particle of rest 
mass to (this assumptions is thoroughly investigated in 
Ref. [3]). The spectral decomposition of Eq. (16) is easily 
obtained from assuming the relativistic mass-shell rela¬ 
tion in flat space, 

M 2 = p 2 + m 2 , (18) 

and by going to momentum space, 

■ 

Mp) = A e 3 / 2^ 3 /4 ’ ( 19 ) 

where p 2 = p ■ p is the square modulus of the spatial 
momentum, and the width A = m p i v jt ~ to. 


For a < 1, one can now write a HWF for each of the 
two horizons in Eq. (15). In fact, the total energy M can 
be expressed in terms of the horizon radii as 


and 



R+ +R- 
2 


R± 


— 


1 ± Vl - a 2 

1 =F a/1 - a 2 


( 20 ) 


( 21 ) 


where M , R + , and i?_ were promoted to operators M, 
R + , and R- related to the corresponding observables. It 
is important to remark that our choice is not unique, be¬ 
cause of the usual ambiguities that emerge when going 
from a classical to the quantum formalism. The unnor¬ 
malized HWFs for R + and R- are obtained by express¬ 
ing p from Eq. (18) in terms of M in Eq. (20), and then 
replacing one of the relations in Eq. (21) into Eq. (19). 
These manipulations yield 


V’H (R±) = A/± 0 (R± — i?min±) 


x exp 


to p R 


2 

± 


2A 2 f2(i±yr 



.( 22 ) 


The step function in the first line above accounts for the 
minimum energy M = to in the spectral decomposition 
of the wave-function (16), which corresponds to 


Rmini = £ P — (l ± v/l - a 2 ) . (23) 

Finally, the normalisations A f± are fixed again by using 
the scalar product in Eq. (7) for both R± 3 . 

The probability density that the particle lies inside its 
own horizons of size r = R± can now be calculated start¬ 
ing from the wave-functions (22) associated with (16) as 


V<±(r < R ± ) = P s (r < R±)Vn(R±) , (24) 


where 


pR± 

Ps(r < R±) = 47t / |^ s (r)| 2 r 2 dr 

J o 


(25) 


is the probability that the particle is inside the sphere 
r = R±, and 


V H (R ± )=4 : TrR 2 ± \MR±)\ 2 (26) 


is the probability density that the sphere r = R± is a 
horizon. Finally, one can integrate (24) over all possible 
values of the horizon radius i?+ to find the probability 
for the particle described by the wave-function (16) to be 
a BH, namely 

nOO 

Pbh+ = / P<+{r < R+)dR+ . (27) 

i„ + 


3 Explicit expressions of J\f± are very cumbersome and not partic¬ 
ularly significant, thus shall be omitted throughout the paper. 
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FIG. 1. Top panel: probability density V<+ in Eq. (24) that 
the particle is inside its outer horizon r = R+ , for m = 2 m p 
(thick lines) and m = 0.5 m p (thin lines) with a = 0.3 (contin¬ 
uous lines), a = 0.8 (dotted lines) and a = 1 (dashed lines). 
Botton panel: probability density P<_ in Eq. (24) that par¬ 
ticle is inside its inner horizon r = R-, for m = 2 m p (thick 
lines) and m = 0.5 m p (thin lines) with a = 0.3 (continuous 
lines), a = 0.8 (dotted lines) and a = 1 (dashed lines). For 
a = 1, the two horizons coincide and P<_ = P<+. 

The analogous quantity for P_, 

poo 

Pbh- = / P<-(r < R_)dR - , (28) 

J Rmi n - 

will instead be viewed as the probability that the particle 
lies further inside its inner horizon. It is already clear 
from these definitions that Pbh- < Pbh+, and it is only 
when Pbh- is significantly close to one that we can say 
that both P_ and R+ are physically realised. 

B. Inner and outer horizon probabilities 

The probabilities Pbh± can only be computed nu¬ 
merically, and we shall therefore display their behaviour 
graphically. 

In the upper panel of Fig. 1 we show the probability 
density that the particle lies inside the outer horizon r = 
R + from Eq. (24) for two values of the Gaussian width 



FIG. 2. Probability Pbh+ in Eq. (27) for the particle to be 
a BH (thick lines) and Pbh- in Eq. (28) for the particle to 
be inside its inner horizon (thin lines) as functions of a for 
m = 2 m p (continuous line), m = m p (dotted line) and m = 
0.5%, (dashed line). For a = 1 the two probabilities merge. 

£ = \ m ~ m _1 (above and below the Planck scale). This 
probability clearly decreases when m decreases below the 
Planck mass, which corresponds to £ increasing above t p . 
A similar analysis is presented in the lower panel of the 
same figure for the probability density that the particle 
lies inside the inner horizon r = P_. It is obvious that 
the probabilities in the second case are much smaller with 
decreasing a. As expected, the two probability densities 
are identical for the extremal case a = 1 (thick and thin 
dashed lines), since the two horizons coincide. 

The probabilities Pbh± are obtained by performing the 
integrations (27) and (28). The plot in Fig. 2 shows these 
probabilities as functions of a for values of the particle 
mass above, equal to and below the Planck mass. When 
analyzing the outer horizon, one notices that Pbh+ stays 
very close to one for mass values larger than the Planck 
scale. However, for m < m p (when the width of the 
Gaussian wave-packet £ > £ p ), Pbh+ clearly decreases as 
the BH specific charge increases to one. Note that this 
probability is not exactly zero even for values of the mass 
smaller than ?n p . For instance, in the case m = 0.5m p , 
corresponding to a width £ = 2 £ p of the Gaussian wave- 
packet, Pbh+ — 0.2 for a considerable range of values 
of cc. It only decreases below 0.1 when a approaches 
one, therefore when the BH becomes maximally charged. 
The situation is very different for the inner horizon. The 
same plot shows that the probability Pbh- starts from 
almost zero for small values of the charge-to-mass ratio 
and increases with a. The larger the mass of the particle, 
the smaller the value of a for which the probability starts 
to become significant. Still, there is a considerable range 
of values of the specific charge for which, while Pbh+ — 1 
thus making the object a BH, the probability for the inner 
horizon to exist is approximately zero. 

Fig. 3 shows the probabilities Pbh± as functions of the 
mass m for a = 0.3, 0.8 and 1. From this plot it be¬ 
comes clear that for smaller values of a, the probability 
Pbh+ starts to increase from zero to one at smaller val- 
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FIG. 3. Probability Pbh+ in Eq. (27) for the particle to be a 
BH (thick lines) and Pbh- in Eq. (28) for the particle to be 
inside its inner horizon (thin lines) as functions of the mass 
for a = 0.3 (continuous line), a = 0.8 (dotted line) and a = 1 
(dashed line). For a = 1 thick and thin dashed lines overlap. 


ues of m. The opposite is true when analyzing Pbh-- 
For the smallest value of the charge-to-mass ratio con¬ 
sidered here, a = 0.3, it is only around a particle mass 
of to ~ 6 to p that both probabilities Pbh+ and Pbh- 
have values close to one, while Pbh+ already increases to 
one around m p . This means that in the range of masses 
between m v and 6 to p , the probability Pbh+ — 1 while 
Pbh- — 0. This mass range increases even more for 
smaller values of the specific charge, but it decreases to 
zero in the maximally charged limit. 

Our main finding is therefore that there exists a consid¬ 
erable parameter space for to (around the Planck scale) 
and a < 1 in which 


Pbh+ — 1 and Pbh- — 0 . (29) 

Whether the particle is a BH or not is dictated by the 
existence of the outer horizon, therefore we interpret 
Pbh+ — 1 as meaning that the particle is (most likely) a 
BH. However, the presence of an inner horizon at r = P__ 
is important in light of the “mass inflation” instability 
and peculiar features of Cauchy horizons. Eq. (29) there¬ 
fore means that the particle is (most likely) a BH, but 
no such peculiarities are expected to occur. 


C. Generalised uncertainty principle 

The uncertainty in the horizon size for a neutral BH 
was already investigated previously. In particular, it was 
shown to grow linearly with the BH mass and lead to a 
generalised uncertainty principle (GUP) in Ref. [2]. We 
can here repeat the same arguments for the outer horizon 
of the charged BH, and obtain similar results. 


We first note that the expectation value 
/»00 

( R +) = 47r / |V> H (R+)| 2 R+dP + 

J Rmin+ 

4(1 + VT^a 1 ) i\ 

2 + e y/neric(l) £ 

= R+(M), (30) 


reproduces exactly the classical expression of P+ in 
Eq. (15) for £ = A m ~ ?n _1 and M = 4?n/[2 + 
e y/n erfc(l)] ~ 1.45 to 4 . From 

nOO 

( R+ ) = 4 7T / |'0h(-R+)| 2 R% dP + 

J Rmin+ 

(l + \/l — a 2 ) 2 (10 + 3 e y/n evic(l)) £ p 
2 (2 + e ypii erfc(l)) £ 2 

~ R 2 + (M) , (31) 

one can then calculate the uncertainty 

AR + = \! ( R+ ) — ( R+ ) 2 csP + ~m, (32) 

which, like in the neutral Schwarzschild case, grows lin¬ 
early with the mass to of the source. This signals the 
fact that the state of such objects would remain QM even 
in an astrophysical regime, where we instead expect the 
horizon has a sharp location, and supports alternative 
models of large BHs, such as the ones in Refs. [6]. 

If we now combine the horizon uncertainty (32) with 
the usual QM uncertainty in the radial size of the source, 


A r 2 = 
~£ 2 , 


|^ s (?’)| 2 r 4 dr — (4'j 


IV’s(r )! 2 r 3 dr 


(33) 


we can finally obtain a total uncertainty 



(34) 


where 7 is a coefficient of order one. The result is plot¬ 
ted in Fig. 4, where it is also compared to the usual 
Heisenberg uncertainty in the size Ar of a state with an 
uncertainty in momentum given by 


/»00 / p OO \ 

Ap 2 = 4 7T / |V’s(p)| 2 p 4 dp — (4 7T / |V>s(p)| 2 p 3 dp 

Jo \ Jo y 


~ a 2 ~ r 


(35) 


We can therefore conclude that the outer horizon behaves 
qualitatively like the neutral Schwarzschild radius. 


4 This mass renormalisation, with M > m, can be easily under¬ 
stood by noting that the source wave-function 'ips contains energy 
contributions from momenta p > 0. 
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FIG. 4. Uncertainty relation (34) (solid line) as a combination 
of the QM uncertainty (dashed line) and the uncertainty in 
horizon radius (dotted line), for 7 = 1. 


IV. CONCLUSIONS AND OUTLOOK 

In this work we transparently applied some basic QM 
methods to the case of the RN BH. In our formalism, the 
location of the horizon is not given by a sharp classical 
value, instead it is described by a quantum wave-function 
with associated uncertainties. In addition, one can define 
a quantity which corresponds to the probability for a 
horizon to be formed. Since the RN BH has the inner 
and outer horizon, the structure of a quantum BH might 
be quite different from the structure of the classical one. 

We first note here that, if one models a BH as a point¬ 
like source (i.e. one very narrow Gaussian wave-function), 
the uncertainty in the location of the outer horizon A R + 
grows with the mass of the source (since £ ~ to” 1 ). This 
is in agreement with many postulated GUPs in the pres¬ 
ence of gravity [ 11 , 12 ], but implies that when the source 
is of astrophysical size, the fluctuations become unac¬ 
ceptably large. This might imply that there must be 
a regime where the GUP gets replaced by the standard 
(non-gravitational) uncertainty principle of QM. In other 
words, the parameter 7 in Eq. (34) becomes very small 
in the semi-classical regime. 

The other (somewhat related) option is that the single 
Gaussian source (16) does not appear to be a sensible 
model for large BHs. The fact that the uncertainty in 
the horizon location does not decrease for increasing rest 
mass to, might imply that a semiclassical behavior can¬ 


not be recovered at all in such a setup [2]. This result was 
remarked in Ref. [5], where it was then shown that mod¬ 
elling a BH as a large number of light constituents [ 6 ] 
does not suffer of this limitation. Because of that, we 
cannot here extend our findings to arbitrarily large BH 
mass straightforwardly. 

We also note that modelling the source with some 
shape different from the Gaussian will not make qualita¬ 
tive differences. For example, using a step function will 
give a result in qualitative agreement with the Gaussian 
of similar width. In fact, one can reliably approximate 
any localised state 4>s with a superposition of Gaussians. 
Given the linearity of the formalism, a superposition of, 
say, N Gaussians of roughly similar mass to will lead to a 
total HWF given by a superposition of the corresponding 
HWFs. In particular, if the width £ of the Gaussians are 
very narrow {£ < £ p ), we can understand what happens 
by simply replacing the superposition with one Gaussian 
of same £ and mass equal to N in. The latter will in¬ 
crease ( Rh ) £, and also increase the probability Pbh- 

However, if £ ~ 1 /V to, the uncertainty A Rh ~ (Rh) 
and the system will never look classical. In contrast, 
if we model the source as N Gaussians of large width 
£^$> £ p , but m <C to p , like in the Bose-Einstein Conden¬ 
sate model of Dvali and Gomez [6], then ( Rh ) ~ Rh 
and ARh ( Rh ) for large N, which does instead look 
classical. If the source has electrical charge, the above 
argument then holds for R + . Since the N constituents 
of the condensate have width £ ~ R + , it immediately 
follows that the inner horizon P_ will have very small 
probability to exist (expect perhaps in a near-extremal 
configuration). 

To conclude, we can speculate that, at least in a quan¬ 
tum regime of BH masses, say to < 10to p , quantum fluc¬ 
tuations around the inner horizon are strong enough to 
prevent the instability expected according to the semi¬ 
classical analysis. More generally, the probability that 
any instability occurs will be as small as the probability 
Pbh- that the source is located inside the inner horizon. 
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